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Abstract: We consider a class of mechanical particle systems interacting with 
thermostats. Particles move freely between collisions with disk-shaped ther- 
mostats arranged periodically on the torus. Upon collision, an energy exchange 
O i' occurs, in which a particle exchanges its tangential component of the velocity for 

a randomly drawn one from the Gaussian distribution with the variance propor- 
tional to the temperature of the thermostat. In the case when all temperatures 
are equal one can write an explicit formula for the stationary distribution. We 
consider the general case and show that there exists a unique absolutely con- 
tinuous stationary distribution. Moreover under rather mild conditions on the 
initial distribution the corresponding Markov dynamics converges to the equi- 
librium with exponential rate. One of the main technical difficulties is related 
to a possible overheating of moving particle. However as we show in the paper 
non-compactness of the particle velocity can be effectively controlled. 

o ■ 



1. Introduction 



Rigorous derivations of macroscopic heat conduction laws from microscopic dy- 
namics of mechanical models require good mixing properties and fast conver- 
gence of initial distributions to the invariant measure(s). For many such systems 
in non-equilibrium, e.g. when a system is coupled to two or more unequal heat 
reservoirs, pure existence of invariant measures is a nontrivial and open ques- 
tion due to non-compactness of the phase spaces. It is relatively easy to envision 
scenarios under which a particle freezes or heats up, which may push initial dis- 
tributions towards zero or infinite energy levels and ultimately violate existence 
of the invariant measures. Tightness arguments are required in order to show 
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that such scenarios occur with zero probability and invariant measures indeed 
exist. Even stronger controls are required to obtain mixing properties. 

Mechanical particle systems coupled to heat reservoirs has seen a renewed 
interest and activity over the last decade. Various examples were introduced in [TJ 
[2HHM71IHJIS] ■ Rigorous results on ergodicity and absolute continuity of invariant 
measures assuming existence for some of these systems have been obtained [IJ 
5,4,1 1]. These results, however, rely on sample paths with infinitesimally low 
probability of occurrence giving no control on times and rates and cannot be 
used to obtain mixing properties or even existence of the invariant measures. 
We present an example of a simple particle system interacting with thermostats 
for which we are able to control more regular sample paths. 

The example is motivated as follows: Consider a system of N non-interacting 
particles at various velocities bouncing elastically off the walls of a bounded 
domain. We assume for visualization purposes that N is very large and the 
system is at temperature To in the following sense: kinetic energies of the particles 
are distributed with the Gibbs distribution with parameter (3q — tjt, i.e. the 
probability that a given particle has kinetic energy dE near E is approximately 
ce-P° E dE. 

Let us introduce a thermostat into the system set at a different tempera- 
ture T\ 7^ Tq such that when a particle collides with the thermostat, an energy 
exchange occurs. That is, upon collision, the thermostat absorbs part of the 
particle's energy, which depends on the angle of the collision, and the particle 
acquires an energy E from the thermostat drawn form Gibbs distribution with 
parameter f3i, where 0i = -^r. Over time, such a system is expected to settle at 
temperature T\ , i.e. the initial Gibbs distribution with parameter /3q is expected 
to converge to the Gibbs distribution with parameter j3\ . The questions of inter- 
est are whether the Gibbs distribution with parameter j3\ is indeed the unique 
invariant measure for the system to which all (or almost all) initial distributions 
converge, and if so, at which rate. 

Now let us add another thermostat at yet a different temperature T 2 ^= T\. 
In this case the system is not in thermal equilibrium. Does an invariant measure 
exist for such a system? Is it unique and if so, do reasonable initial distributions 
converge to it and what are the rates of convergence? Similar questions may be 
asked in the presence of more than two thermostats at different temperatures. 

In the absence of particle interactions the system with many particles is 
simply the product of one particle systems. The dynamics is described by a 
continuous-time Markov Process, which is deterministic apart from collisions 
with thermostats and upon a collision of a particle with a thermostat, a random 
perturbation occurs. This degeneracy of the Markov process allows to restrict 
the study to the discrete time dynamics on the collision manifold. 

For the resulting discrete-time Markov chain we show that, under certain geo- 
metric assumptions, there exists an invariant measure and it is unique (ergodic), 
absolutely continuous with respect to Lebesgue measure, and mixing with ex- 
ponential rates. We also conclude that reasonable initial distributions converge 
to the invariant measure exponentially fast with control on the rates. It follows 
that for the original Markov process there exists an invariant measure and this 
measure is unique (ergodic) and absolutely continuous. Mixing and convergence 
of initial distributions to the invariant measure for the Markov process do not 
follow directly since under some scenarios particles may move extremely slow. 
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Though ergodicity guarantees that they will eventually speed up, bounds on the 
times must be obtained in order to show mixing. We leave the investigation of 
mixing properties of the Markov process for future work. 

The proof of mixing properties in the discrete case uses general state Markov 
chain machinery, in particular, Harris' Ergodic Theorem [51110] . The theorem 
requires two things: to produce a non-negative function V on the phase space 
which, on average, decreases geometrically under the push forwards of the dy- 
namics, and, given such a V, to show minorization or Doeblin's condition on 
certain level set of V. The first condition guarantees that the dynamics enters 
the 'center' of the phase space, a certain level set of V, with good control on 
the rates; and, once at the 'center', coupling is guaranteed by the minorization 
condition. Those two conditions imply existence and uniqueness of the invariant 
measure with exponential mixing rates and exponential convergence of reason- 
able initial distributions to that invariant measure. 

We believe that most of our results can be generalized to the case of ther- 
mostats of a general smooth convex shapes rather than disks. It is also interesting 
and tempting to extend the methods of this paper to the 3-dimensional case. We 
are planning to address these and other open questions in the near future. 

We describe precise settings in section^ then state the results and outline 
the proofs in section [3] In section H] we describe the potential V, and in section 
[5] we establish the minorization condition. 



2. Settings 

Let r = T 2 \\jP =1 Di be a bounded horizon billiard table with finitely many non- 
intersecting circular obstacles D±, • • • , D p having radii Ri, ■ ■ ■ , R p . An example 
of such configuration is shown in Fig. [TJ Particles move around in r. Apart 
from collisions with the boundary <9P, the particles move freely with constant 
velocities and do not interact with each other. Each obstacle Di, 1 < i < p, 
is assigned certain parameter /3i which represents the inverse temperature and 
plays a role of a thermostat. 

Upon a collision of a particle with a thermostat Di, let v~ — be 
the decomposition of the particle's velocity into the normal and the tangential 
components with respect to the boundary of the thermostat ODi. After the col- 
lision, the normal component of the velocity changes sign, = — uj, while the 
tangential component is absorbed by the thermostat and a new tangential com- 
ponent vf is drawn from the distribution \J~^-e~P iVt dvf . The outgoing velocity 
has a decomposition v + = {v\,vf). This type of energy exchange was used in 

m. 

Since particles do not interact with each other, we can reduce our discussion 
to studying the system of only one particle. The Phase Space of such a systems 
is 

Q = {{x,v) : x G r,v G R 2 }/ ~ 

where ~ is an identification of points on the collision manifold that correspond 
to incidence and reflection with the choice of the old, incoming, velocities v~ = 

(«!»«*")• 
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Fig. 2. Discrete Dynamics Coordinates 
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Then the dynamics is described by a Markov Process <P T on fi, which is deter- 
ministic between collisions with thermostats with random kicks at the moments 
of collisions with dr. 

As usual for such systems, it is convenient to pass to a discrete dynamics on 
the boundary dT of r. Note that at the moment of collision we chose to keep old, 
incoming velocities. Given such a choice, to describe the discrete dynamics it is 
enough to keep track of the sequence of collision points r n G dT (parameterized 
by arc- length) and the absolute value of the normal velocity v± (n) at the moment 
of collision. Using randomly generated tangential velocity v t (n) drawn from the 

probability distribution \j\eT^ dv t , one determines the outward velocity and, 

hence, the next point of collision r n+ \ and the next normal velocity v±(n + 1). 
This procedure defines a Markov chain <P on Q = dr x [0, oo). 

For our future analysis it would be more convenient to use a geometric de- 
scription similar to one used in billiards. Namely, instead of the random vari- 
able v t (n), we shall consider the random variable ip n , an angle relative to the 
normal vector at the collision point at which a particle leaves after collision 
n. It is an easy calculation to show that ip n is drawn from the distribution 

p v±(n )(<P n )d<p n = V^SS^ e ~ ft " i( " )tan2(v "W™, which depends on ft asso- 
ciated with dDi 9 r n and v± (n) . We shall also use notation ip' n for the incoming 
angle at the point of collision r n +i. Obviously, r n +i, vj_(n + 1), and ip' n are com- 
pletely determined by r n , v±(n), and ip n . Note that v±(n + 1) = 

Denote the transition probability kernel of $ by V((r, v±), •), i.e. V((r, v±), A) = 
P(<£„ € A\<P n -i — (r, v±)). We will also use the notation V* and V* for the op- 
erators defined on both the set of bounded measurable function and the set of 
measures of finite mass by 

(Vf)(x,v ± )= [ /(/,«l)P((r,v L ),d(r , y_ L )) 
Jn 

(V^)(A)= f V((r,v±),A)n(d(r,v ± )) 

Note that the transition probabilities are degenerate: V((r, v±), •) = V*6( rtV± ) is 
supported on a family of one dimensional curves in the two dimensional phase 
space fi. 

Our interest lies in investigating the questions of existence, uniqueness, abso- 
lute continuity w.r.t. Lebesgue measure m, ergodicity, and mixing properties of 
the stationary (invariant) measures for the Markov Chain # and the associated 
Markov process 3 T . 



3. Results 

We show existence of the invariant measures in the non-equilibrium situation 
simultaneously with the geometric ergodicity for the Markov chain <P. 
The main result is given by the following Theorem: 
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Theorem 1. The Markov Chain <P admits a unique absolutely continuous in- 
variant probability measure fi. Furthermore, there exist a non-negative function 
V on Q, as well as constants C > and 7 <E (0,1) such that 

sup \V n ((r,v 1 _),A)-^(A)\<Cra + V(r,v ± )) 
ago 

Corollary 1. The invariant measure of the Markov chain <P is exponentially 
mixing, i.e. for any Borel A, B G Q, 



sup I / V n ((r,v±),A)dn- n(B)n{A)\ < Cf l / (1 + V(r, v±))dfx < De~ an , 

AGO J J 
B B 

for a — — 111(7) an d D = C J(l + V(r, v±))d[i < 00. 

n 

Corollary Q] follows from Theorem [1] once J(l + V(r,v±))dfj, < 00 is estab- 

n 

lished. This is one of the results of the Theorem 3.2 in [6]. 

Corollary 2. If the probability measure v on Q satisfies J n {^ + V(r,v±))di> < 
00, then 



\\P>-»\\ <CT / (1+V(r,v ± ))dp, 

JQ 

where || • || is a bounded variation norm. 

By constructing a suspension flow over <P we obtain existence, uniqueness and 
absolute continuity of the invariant measure for the Markov process both for 
one and for many particle systems. 

Corollary 3. There exists a unique (ergodic) absolutely continuous invariant 
probability measure for the Markov process <P T 

Remark 1. In the equilibrium case, i.e. when /3i = fti = ■ ■ ■ = fl p = ft the 
invariant measures can be written down explicitly. 

Lemma 1. The measure ji with density 



2/? 2 

\dr\ 



is invariant for the Markov chain <P. 



We will prove Lemma [T] at the end of section [SJ 
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3.1. Idea of Proof of Theorem\J\ To prove Theorem [T] we will use the following 
formulation of the Harris ergodic theorem for Markov chains also referred to as 
Geometric Ergodicity Theorem. 

Theorem 2. Jfjjj Assume 

Potential Condition . 

There exists a function V : fl — >■ [0, oo), K > and 7 G (0, 1) such that 

V*V(r, v±) < jV(r, vj_) + K 

for all ir,V±) G fl and 
Minorization Condition . 

There exists a probability measure v supported on C, N and rjN G (0, 1) such 
that 

inf V N ((r,v ± ),-) > VN v(-), 

where C = {(r, v±) G fl : V(r, v±) < S} for some S > 2K/(1 — 7) where K 
and 7 are the constants from the Potential Condition. In addition, to ensure 
aperiodicity, we require that the same holds for N+l and some rj^+i G (0, 1). 

Then <P admits a unique invariant measure fx. Furthermore, there exist C > 
and 7 G (0,1) such that for all (r, G fl 

II (V* n f)(r, wx) - M(/)(r, «x)|| < C^Wf - 

for all f such that \\f\\ < 00, where \\f\\ = sup i^fe^X) 

In particular, for any A C fl, let / = | q' ^ ^ , then 

|^"((r, Vx ),^)-M^)l < sup \r* n (f) - rtf)\ < cy 



so that 



V(r,w±) + 1 " (ri „ x) y(r,u_L) + l 



sup |P n ((r, » x ), A) - < C 7 n (l + T/(r, uj.)), 

Acn 



which is exactly a conclusion of Theorem Q] We chose a weaker formulation of 
the theorem in order to make the presentation more intuitive. 

To prove Theorem Q] we need to show that the Markov Chain <P satisfies the 
conditions of the Harris Ergodic Theorem and that the invariant measure we 
obtain is absolutely continuous with respect to the Lebesgue measure on fl. We 
construct the potential V in section 0] and prove the minorization condition on 
a level set of V and absolute continuity in section [SJ 
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4. Potential V 

Proposition 1. There exists a junction V : fl — > [0, oo) and constants K > 
and 7 G (0, 1) such that 



V*V(r,v ± )= / V(r , ,v , ± )p Vl (<p]Pi)dtp<>rV(r,v ± ) + K, V(r,vx)Gtt, (1) 



and r' and u' are £/ie position and the normal velocity at the next collision given 
r, v±, and ip. 

The density p vx (ip; fy) implicitly depends on r through the inverse temper- 
ature of the corresponding thermostat. Since we have only a finite number of 
thermostats, the ft are bounded above and below and do not play an important 
role in the asymptotic analysis. To simplify the notation we will write p v± (tp) 
instead of p v± (tp; ft). 

4-1- Heuristics. The inequality (TTJ) ensures that the values of V at the random 
images of (r, v±) are, on average, smaller than V(r, vj_) when V(r, v±) is large, 
and, on average, smaller than a constant when V(r, v±) is small. This implies the 
dynamics enters the 'center' of the phase space, represented by some level set of 
V, regularly with tight control on the length of excursions from the 'center' 
110] . Since we expect that particles reach very low or very high velocities with very 
low probabilities and would like the system to satisfy the minorization condition 
on the center, a natural physically meaningful candidate for the center would be 
a set of states with moderate v± velocities, e.g. C = {(r,v±) : vf m < v± < i;™ ax } 
for some < vf m < vf aK < oo. Thus we are looking for a potential V that takes 
0(1) values in C and blows up as v± — ¥ and as «i — ^ oo. This is essential in 
order to satisfy the minorization condition. 

First let us investigate heuristically the mechanism that ensures that a po- 
tential of this kind should, in principle, satisfy the inequality in Prop. 14.11 As- 
sume for simplicity that V depends on v± only and C is a level set of V, i.e. 
C = {(r,v x ) : vf u < v± < vf ax } = {(r,v x ) : V(v±) < A} for some A. We 
are interested in the change of potential we can expect if we start with low, 
moderate and high v± and iterate one step forward along some likely random 
trajectory. The velocity of the particle between collisions with thermostats is 
v = vj_/cos(ip), where ip is the angle drawn from p v± (ip)dip, and the normal 
velocity upon collision is v'^ — wcos(<//), where tp' is the angle of incidence of 
the particle that originated at (r, tp) . 

— If v± is very small, with large probability vt is the main contribution to the 
overall velocity v and v 1 ^ k, cos(tp')vt- Again, with large probability, tp' is 
bounded away from zb^ so that v'j_ is of moderate range. Thus with large 
probability the value of V drops, i.e. V(v'j_) < V(v±). 



tt/2 



tt/2 



where 
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— If vj_ is 0(1), with large probability the vt contribution is 'comparable' to v±, 
so v is of the same order of magnitude as v± . Again, with large probability 
ip and ip' are bounded away from and thus v'j_ is of the same order of 
magnitude as v± and V(v' ± ) = 0(1). 

— If v± is very large: with large probability v t is negligible and p> ~ 0. Thus 
v' ± ss cos((p')v± < v± provided cos(p') is bounded away from i.e. with large 
probability V(v' ± ) < V(v±). 

We conclude that, in principle, a potential V taking 0(1) values in C and 
tending to infinity as v± — > and as v± — > oo should satisfy the inequality in 
Prop. 14.11 The escapes to infinity should be fast enough to ensure large enough 
drops of averaged value, though not too fast to ensure integrability of the left 
hand side. There are also several potentially dangerous geometric locations: when 
v± is very small near locations r G dr such that <p — — p' — ±-| and when v± 
is very large near r € dr such that ip = (p' = 0. The heuristic large probability 
arguments fail for those. In the next subsection we will present a potential V 
that depends on uj_ and show that it satisfies the inequality of Proposition for 
all small, moderate v±, and large v±, even in the aforementioned dangerous 
geometric locations. 



4-- 2. Construction of the potential V . A reasonable guess for V(v±) is 

V(v±) = { vl a , v x < v'l in ; (2) 
A, vf'" 1 <v±< vf 3 *. 

where e < /3 m i n = min{/3i, • • • ,/3 p } and < a < 2. Let us fix a = — 1. We 

will fix e later. We will also assume that e e ( v± ' = (w" lln ) _1 = A so that V is 
continuous. 

We are interested in establishing that for some 7 > 1 and K < 00 



< [ 



V(v' ± )p v± (<p)d<p (3) 



^^]-V„(^)^+ / exp(e[^^i] 2 )^M# +J 4 / p.A<P)*P 



CO$((f) ^ J cos(^) 

< jV(v ± ) + K 

4-2.1. Controlling small v±. In this subsection we are going to check the in- 
equality for a potential @ given (r,V±), where vj_ < u™ ln for some v™ m small 
enough. When the random map ^ is applied, we first draw ip from the distribu- 
tion with density p v± (p) and then jump to (r', p' , v' 1 ) 1 where (r', ip') is the image 

of (r, (p) under the billiard map and v'^ = 7^^y u _L- As ip sweeps from — ^ to 
the billiard map image (r',(p') sweeps through several thermostats (we only 
count those with positive measure of (r', <//)-images). These thermostats can be 
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divided into two groups: 'outer' thermostats, roughly speaking those that in- 
tersect two rays originating at (r, ±5) and 'inner' thermostats, all the rest. Of 
course, some thermostats in this description may repeat: nevertheless we may 
classify the thermostat differently depending at which angle we are looking at 
it. 

Figure |3] gives an idea of how ™ h s ^^ behaves as tp sweeps from — ^ to ^ in 
a typical case. Each curve represents a thermostat. The figure does not specify 
which thermostat blocks another, so there is an overlay of the values and we 
will be integrating over all the overlays since we do not make any assumptions 
on geometry; of course, only one branch works for each value of tp in the real 
setting. 

Fix < 7 S < 1. We would like to split the first integral in ([3]) into two parts 
depending whether cos ( y / < — or cos ( y ) > J_. 

' o rosier! -v. coslco) — 



cos((p) j s cos(ip) 

r/2 

r COs(^>-L!-l. I Pi v ±- „-A itan 2 («p) 



COs(p) V 7T cos 2 (tp) 



dp 



-ir/2 



COS(p) V 7T COS* (tp) 



cos(p') x f% Vj_ 



cos(^) 
COs(y) V 7T cos 2 (p) 




-PiV 2 ± tan 2 ((/3) 



dp 



7a«x +V / 2TT 2 — 4 

1 ■ ■- ' COS^(p) COS^((/?) 



where /3 max = max{/3i, ■ • • , /3 P }. 

To estimate the second integral in (HJ, we need to express cos(p') in terms of 
tp for each thermostat in range. Given a thermostat Di let </?o be such that corre- 
sponding p' = — § if we assume that all other thermostats are invisible, i.e. the 
particle can pass through them. See Fig. [U For some of the 'outer' thermostats 
such po will not exist: in this case we have to choose tpo with corresponding 
tp' = ^ and perform a mirror-image parametrization of cos(p'). The computa- 
tion is essentially the same and thus omitted. 

By sine theorem (see FigH]) 

R R — To tan(p — po) 

cos(</5 — tpo) sin(— tp') 

Thus 



sin(— tp') — cos((y5 — po) — ^ sin(p — po) = cos(y>) cos(<^o) + sin(y) sin(</>o) 

it 

-— sin(<p) cos ^o) + cos ^) sm((p = ^= ^— 

R R Vw 2 + ly/wl + 1 



Ergodic properties of random billiards driven by thermostats. 

3.0 h 



11 




1.5 



Fig. 3. cos (' p j for all visible scatterers 




Fig. 4. Parametrization of ip 1 

here we changed variables w = tan(t^) with w$ = tan(^o)- 

cos 2 {if 1 ) _ (w 2 + l)(w 2 + !)-(! + wwg - f w + f w )(l + w^p - g w + %wp) 

COS 2 (ip) Wq + 1 



2^ - ^* 



« 2 W 2 



w 



R R 2 w Q-t- R W 



w - 



wl 



aw 



2bw 



For 'inner' thermostats, which we define as thermostats with corresponding 
a < 0, the above expression is a concave down parabola. We are interested in 
the values between the roots, where the smaller of the roots is wo. The maxi- 
mum value at each parabola is at least 1 since for some ip, the particle hits the 



thermostat radially, i.e. with cos(tyj') = 1, so that 
the maximal value of the parabola if aw 2 



2 (<p') _ 



— i — — — irr > 1- Let I be 
2bw + c = 4y has no real solutions 
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and I = \ otherwise. Denote the smaller or the only root of aw 2 — 2bw + c= I 

Is 

by wi. Then 

l(w — w ) o 1 r n 

- < aw — 2bw + c< — on [w ,wi\ ana 

{wi - w ) 



Wl Wl 



[ (aw 2 - 2bw + c)- x ' 2 dw < I' ( 1 } W Wo \ )-V 2 dw = 2(wi - w )rl 
J J (wi - w ) 

Wo w 

It is left to estimate (wi — wo). Simple computation using the values of coef- 
ficients a, b, and c ensures that Vb 2 — ac — jt and that 



(wi - w ) 



R 2 + °' li 



T 

has no extremum points for a e [— j^j- , 0] . Thus 

(wi — w ) < max{ — , lim(wi — w )} = — 

T a^O r 



We conclude that 



f (aw 2 ~2bw + c)- 1 / 2 dw<2— I 

Tn 



Rl,-L . 2i? max 

2 < . 



Ml ' J I Li I I ~ 

Wo 



where i? max = max{i?i, • • • , R p } and T m ; n is the minimum distance of flight. 

The situation is clearly symmetric for another set of roots of aw 2 — 2bw+c = 
and aw 2 — 2bw + c = I. 

For the 'outer' thermostats, when a > 0, aw 2 —2bw+c is a concave up parabola 
or a line with slope 2^ (case a = 0) and we are interested in integrating values 
with w > w - The derivative of aw 2 — 2bw + c at w is 2^ and 2^(w — w ) < 
aw 2 — 2bw + c on [wo, oo]. Let wi = w n + 9 . Then aw 2 — 2bw + c > \ on 
(wi , oo] and 



Wl Wl 

Rn 



f (aw 2 - 2bw + c)~ 1/2 dw < [ (2^(w - w ))- 1/2 dw = — <- 

J J R T -f s T n 

w a w a 

Let T max be the maximum distance of flight on the bounded horizon table. 
Then for each point r e dT the number of 'visible' thermostats cannot exceed 
the maximal number S of thermostats in any disk of radius r max if we view 
the positions of thermostats as a periodic configuration on the plane. Therefore 

cas Jf\ is less than \ on at most 25* intervals and the integral f cos nf ) dw on 

each is < 2fimax . Therefore 

— T min 7 s 
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t/2 



-it/2 



! cos((p')v± . 
cos(<^) ' 



-fiiv± tan 2 (v?) 



< IsV 7 1 + 



A, 



r cos 2 (<//), i dyj 



COS 2 ^) COS 2 ((/?)' 



, —1 . 2i? max / Pijitui L -, 
< 7 s w ± H a / 2S. 



TminTs » 7T 



(5) 



Note that estimate ([5J works equally well for all values of v±; we will use it 
later for large Uj_ estimates. 

The second integral in the equation (|3]) can be bounded by using the fact that 
K) 2 <vl+v 2 t : 



cosjyyL 2 



dvt 



< e e 



- e)v *dv t 



(6) 



Thus 



' C0S(<y9j ' 



cos(w')«j, 2 . . 
r^H / p v± {ip)dip 

COS(ip) ' 



^ —1 2i? max /pmaji -, 

< 7 S w ± H a / 25- 



T"min7s V 7T y ft - e 

for any wj_ < u™ ln for some w™ ln and if s < oo 



exp[e(uf n ) 2 ] + i4 = 7^7/ + X s 
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J^.2.2. Controlling large v±. To deal with large v± we first note that for any 
« -d. 



V(v' x )p v± (<p)d<p 



< 



Ox) 1 Pv ± (f)dip+ / e £( ^ )2 p v± (p)dp + A / p v± (<fW 



Tminls V 7T J V 7T 



<se tv± + e e(v ^> J—e-^dvt, (7) 

for any > ii™ ax , where ^™ ax = ?;™ ax (<^) is large enough. 

Given e > the estimate © gives a bound by J ^i^e £ ' u "" x when < t>™ ax . 

However, this estimate is too rough to produce uniform 7/ < 1 for large v±_ > 
v™, so we need to be more careful here. 

For any <5 > 0, and A large enough, using Chebyshev's inequality, we have 



\v t \>A 



TT 



Pi-e J V TT " V A - e 2(A - e)A 2 < 2 6 ± ' 

Note that A can be chosen independent of e provided e < eo < min{/3i, • • • , /3 P }. 
In the remaining the integral 

e «f^K e -^\ dvt = J e<^p v ^)d V , 

when v± > v™ ax with -y™ ax large enough, we only integrate over values of ip ~ 0. 
So our new perpendicular velocity v' ± = c °^^ v±_ ~ cos(y>o)i;_L depends for 
the most part on the landing angle ip' which corresponds to tp — 0. For most 
locations r € dr, cos(<y9p) is uniformly bounded away from 1. 

Before proceeding with further estimates, we are going to identify trouble- 
some geometric locations. For this purpose we are going to lift the periodic 
configuration of thermostats from a torus to a plane. 

Let / be a line segment connecting the centers of two thermostats Di and 
Dj such that I does not intersect with any other thermostats. Then if a particle 
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originates at r = I fl dDi in the normal direction, i.e. with ip = 0, it must hit 
Dj at r' also in the normal direction, i.e. with ip' = 0. In addition, the normal 
velocity does not change along such a trajectory, i.e. v'j_ = v±. Suppose now that 
initial v± is extremely large. Then the tangential velocity vt acquired before the 
flight is, on average, tiny compared to v± and the particle lands very close to r' 
with <p' w 0; i.e. v'j_ w v± and we do not get much of a drop in v±. And it is clear 
that the higher the initial velocity v±, the smaller is the drop. This indicates 
that locations like r and r' are 'troublesome' and need to be treated with extra 
care. 

Given a thermostat D i7 consider all line segments Ij on the lift to the plane 
connecting its center to the centers of other thermostats, possibly intersect- 
ing other thermostats, and of length at most 2i? max + r max , where i? max = 
max{i?i, • • • ,R P } and r max is the maximal length of collision-free path. Note 
that there are only finitely many such segments Ij since there is an upper bound 
on the length. Let rj = Ij n dDi. If a particle originates from r £ dDi in the 
normal direction denote the point on next collision with a thermostat by r' 
(tangential collisions count as collisions). 

Let Ari be small enough such that for any j the intervals [rj — Ari,rj + Art] 
do not intersect and for all r £ [rj — Ari,rj + Art], all r' belong to the same 
thermostat unless Ij touches a thermostat tangentially. In the latter case we 
require the same for all r £ [rj — Ar i7 rj) and for all r £ (rj,rj + Ari\. 

Let Ar < mim{ R ™™ Ari } and denote by TZ the union of all r/s over all 
thermostats Dj. In addition, consider all segments / connecting the centers of 
thermostats but not intersecting any other thermostats except tangentially. De- 
note the collection of possible intersections of such / with dT by {ru}- Note that 
{rfc} C TZ and that if r £ (r^ — Ar, + Ar), for all <p, such that \ip\ < R r , r' 
belongs to the same thermostat as r' (since (r^ + 2(r — rk))' does). 

Lemma 2. Suppose 1-segment path originates from (r,v±). Then we can choose 
Ar < Ar small enough such that 

— ifr& (rfc — Ar,rk + Ar) for all k and \ip\ < ^ r , then \<p'\ > 2 ^ r . 

— if r £ (r/j — Ar,rk + Ar) for some k and \ip\ < ^~ rfc ^ ; then \ip'\ > ^ rfc 

Proof. Choose Ar < Ar such that sin( ^ r ) > 2 Ar . We are going to focus on 
the first statement. The proof of the second statement is essentially the same 
noting that the situation with tangential collision would only give a better esti- 
mate. 

Denote the line segment from (r, ip) to (r', p') by 7, the thermostat containing 
r by D, the thermostat containing r' by D', radius of D by R, and radius of D' 
by R' . Let I be the line segment connecting the center O of the thermostat D to 
the center O' of the thermostat D' and let rj — I n dD and r\ = I n dD'. Note 
that |/| < 2i? max + r max since |— y j < r max . By assumptions on Ar, we conclude 
that \r — rj\ > Ar. 

Let J be the line segment parallel to / originating at r and ending at some 
r'| . Denote the angle J forms with the normal to dD by p and with the normal 

to the disk of next collision by p' . Then |<^o| = - fl ' J > ^ > |^|- Therefore 
r' £ dD' and \p'\ > \p' \. 
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Then > \ 9 ' \ = {Lr'fi'r'A > | sm^O'rf,)! = |«S£^IiM| > ^ > 



2i{ m „ ■ 

.max \ AR 



If r ^ (rfc — Z\r, rfe + Z\r), then for v t < A and v±_ > vj ax > ^™ ax large 

v ± - R n 



enough, tp < tan((p) < — < J^ r and 



«) 2 = (vl + v?) cosV) < «i - sin 2 (-^-),i + A 2 < vl - avl, 
where a = j sin 2 ( 2J ^ r ). 



Given 5 > and /l satisfying ^/ -g%^ 2 (/3 -e)/i 2 < f we 



V 7T 



7T J V 7T 

Kl>/i M</i 

< 5 -e tv ^ + e~ eav ^e e < <d + l- 8)e ev ± = (1 - 5 -)e tv ^ 

provided that v± is large enough so that e~ €av± < 1 — 6. 

Therefore we need to focus our attention on studying the dynamics in regions 
(r-fc — Ar, rk + Ar). Let r £ (r^ — Ar, r^ + Ar) and r' be the image of r given the 
tangential velocity vt < A. First let us note that we can ignore the situation when 
intervals / connecting centers of thermostats intersect some other thermostats 
tangentially since it would only give us a better estimate. Then if v± > y™ ax > 
A ^" ax all appropriate images r' lie on the thermostat D' k . 

Let 

rv± , r'v±_ 

y= ~R~ and v = ~W 

Lemma 3. Given any A > and A > 0, there exist Y = Y(A,A) and v^ x = 
v^ x (A,A) such that 

W^f-vlK-A 
provided that \y\ > Y , \v±\ > ti™ ax and \v t \ < A. 

Proof. Pick Y = max{R max A, 8R niax (A + A 2 )}. Then if \y\ > Y, \r - r k \ > £. 
And by Lemma[H since \(p\ < \ tan(ip)\ < -A- < R ^ v , 

W) 2 ~vl<- sin 2 ( — ^—)v\ + /l 2 < + A 2 <-A 

provided that v±_ > i>™ ax is large enough. 
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Then if we choose A such that e eA < 1 — 5, for \y\ > Y and \v±\ > w" lax 



OO 

/ 



'Pie-^dvt < £e CT * + e~ f - A e f -< < d + 1 - o)e e ^ = (1 - ^)e e 
7r 2 2 2 



It remains to treat \y\ < Y case. 
Lemma 4. Suppose \v t \ < A and \y\ < Y. Then 



-{R' k y R' k s 4 L fc J 

+ [2dfi fe + fi^ + rf 2 + (i?;) 2 + fifefife]^ +£(«±), 

w/iere £ — > as uj_ — > oo uniformly in {(y,v t ) : < F, |«t| < yl}. fiere fife 
and fi fc are radii of D k and D' k respectively and d — dk,k' is the distance between 
Dfc and D' k . 



Proof. From geometry 



fife + d d , 

r H «t + (Jyr ) as r — > 



fifc ' «_L 

Then, since |y| < Y 

, fifc+d d 1 
y = — ^7 — y + -^j-v t + 0( — ) as vj_ Y oo 
^fe ^fe U -L 

Note that both 0(r 2 ) and 0(777;) are uniform in {(y,v t ) : \y\ < Y, \v t \ < A}. 

Let us orient the appropriate line segment I connecting centers of D k and D' k 
vertically and let 9 = and 9' = j^- be the angles r and r' forms with /. Then 

Vvert = v± cos(6») - v t sin(6>) 
Vhor = v± sin(0) + v t cos(6>) 

and 

v' ± = v vert cos(O) - v hor sin(6>) 
= v± cos(6>) cos(0') - v t sin(0) cos(6>') - v± sin(0) sin(0') - v t cos(8) sm(9') 

= «x(l -\^T~ \^T~) ~ vt — 1 ~ v±(——) - vXl + 0{\) 
2v z L 2 vj_ v± v±v± v± vj_ 

, 1 (y 2 + W? + 2v t y + 2yy' + 2v t y') 1 

= 1>± + h 

2 v± vj_ 

So 

«) 2 - v\ = (y 2 + {y' f + 2v t y + 2yy' + 2v t y') + O(^) 

fi fc + d d 2 R k +d d d 2 2d 2 

H k H k H k H k K H k) H k 
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| 2dv t + 2(R k + d)dv t 
R' 



(K) 2 



y + 2v t y + 2v t y 



R k + d 
Rl 



[R k + R' k +d] 



2 V 



(K) 2 



£{v A 



2d 



— }v 2 + [2dR' k v t + 2{R k + d)dv t + 2v t {R' k f + 2v t {R k + d)R[ 



R 



+ [R k + R' k + d] 2 j^ + £(v±) = + ^-}v 2 + [R k +R' k + d} 2 ^ 

+[2dR' k + R k d + d 2 + (R' k ) 2 + R k R' k ]^j + £(v±), 
where £(v±) = 0(v±) as v± — > oo uniformly in {(y,v t ) : \y\ < Y, \vt\ < A}. 



We are now ready to estimate e ev± J e^'"^ 2 y— e P* v * dvt for y < Y and 

\i>t\<A 



exp(-Pivt +<v' ± Y -ev±) 



IPi , ]tt , d 2 2d 2 



[(R' k + df + dR k + R k R' k ]^L e - [R k + R' k + dfjj^e + £(v ± )) 



yRe , 



§V§«p(-fl(«*-. 



if) 



c M-T^[{Rk + K + d) 2 e + 8(v x ) 



R 2 



(R' k ) 2 (3. 



where ft = ft(e) = (3 t + + e# and R = {R' k + d) 2 + dR k + R k R' k . 

\ k ' k 

Let e < eo be small enough and {j™ ax large enough such that 



[R k + R' k + d] 2 e + £{ Vl _) 



R 2 



(Km 



e 2 >0 



for v± > v™ 3 -* (can do this by Lemma HJ. Then 



l«t|<-A 
Let 



7 



max 



/A 



max 



A 



{Rfc.R'fcA^.ft} V A {«*.-RL.«ifc,fc'./3<} V A + 



- 2d 



Note that 7 < 1. Choose S and <r so small that 7 + (5 + q < 1 and 1 — | + <; < 1 
and then choose A so large that 



ln 1 _ 6 ) yl 2 < f . Let 7/ = max{7 
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5 + q, 1 — | + <;}. Choose v™ ayi > {j™ ax that satisfies Lemma [3J and all the other 
assumptions we made. Then for v± > 



V(v' ± )p v± (ip)d(p < 7 iexp e ^ = tjV(«±). 



Let 7 = max{7 s ,7i} and K — max{i\~ s , J ^^■" i ^^ e e ( t, T ax ) 2 |. Then 



tt/2 

P*V r (»±)= / ^KK(^<7^M + ^ V(r,« ± )el2. 

-7T/2 

5. Minorization Condition 

For any u™ ln < u™ ax , let C = {(r, : v™ m < v± < u" lax } and let v be uniform 
probability measure on C. 

Proposition 2. There exist an integer N > and t]n > such that 

inf V N ((r,v ± ),-)> m u(.) 

(r,v ± )EC 

In addition, there exists tjn+i > such that the same holds for (N + l) st push 
forward. 

We are going to prove Prop. [2] in three steps: (i) Prop G2 which guarantees 
existence of a 'regular' iV-step sample path from any (r, v±) € C to any (r', v' ± ) € 
C, where C = {(r,w_L) : v™ in < v±_ < w^ ax } for some 6™ in < v™ in and w^ ax > 
w™; (ii) Prop. |H which gives bounds on the densities when pushed forward 
along regular paths; and (iii) Prop. 02 which guarantees that for any (r,v±) € C 
we can 'acquire' density in two steps and the resulting location (r", t>" ) € C. So, 
originating at (r,v±) G C, we acquire density around (r", ) € C by Prop [5] and 
then push it forward along the sample path from Prop [3J and Prop. U ensures 
that we get a uniform lower bound on the density at the final point (r', 

To state Prop. [31 [H and[5j we need the following definitions. 

Definition 1. A projected (particle) path is a continuous curve 7 : [0, 1] — > r, 
s 1— > 7(s) i/iai consists of a finite sequence of straight segments meeting at dT 
with 7(0), 7(1) G dr. 

Remark 2. Note that a projected path is allowed to have any 'reflections' off 
the boundaries of the thermostats UjL^dDj. An example of a projected path is 
shown in Fig[T] 
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If we start a projected path with certain normal velocity v±, then all subse- 
quent normal velocities upon collisions are completely determined by the path. 
More precisely, at k th collision 

yk = cos(y? / 1 )---cos(y fc ) ^ 
^ cos(^i) • • • COS(^fc) 

where ipi , • • • , ifk and (p[, ■ ■ ■ ,(p' k are the angles of reflections and incidences for 
the path. Note also that multiplying the initial v± by a certain factor multiplies 
all subsequent v\ by the same factor. 

Definition 2. Let 7 be a projected path that starts at r and ends at r' . A sample 
path a given 7 and v± is a path in the phase space that starts at (r, v± ) with po- 
sitions updated according to positions of the path 7 and normal velocities updated 
as above. 

Note that given 7 and v± , a is unique and to each sample path a corresponds 
a unique projected path. 

Definition 3. Given S > 0, Ar > 0, and Atp > a projected path 7 is (S, Ar, Aip)- 
regular if 

— All the incoming and outgoing angles are bounded away from ±-| by S, i.e. 
all(p,<p' e(-f + £,§ -8). 

— For each segment of 7 from (rk,(fik) to (rk+i,<p' k ), any segment originating 
fromfk G (?*& — Ar, r^ + Ar) with angle ipk G (Vk — Atp, (fik + Aip) has incoming 
position fk+i on the same thermostat as rk+\ and angle <p' k G (— -§ + 5, § — S); 
in addition, the same holds for each segment in reverse direction, i.e. from 
( r k+i,p' k ) to ( r k,<Pk)- 

Proposition 3. There exist N, So, Aro, and Atpo such that given any (r,v±) 
and (r',v' ± ) with u™ m < v±,v' ± < v™ ax , i.e. (r, v±), (r' , v' ± ) G C, there exists a 
sample path from (r, v±) to (r', v'j_) such that the corresponding projected path 

1. is (So, Arg, Aipo) -regular. 

2. has exactly N collisions. 

Here vf in = cos 2 (f - S )vf in and C? ax = "/T , » . 

_L \2 u/ -L -L cos^ ( £ — do ) 

We are going to prove Prop. [3] in three steps. The first two are covered by the 
following Lemmas. 

Lemma 5. There exist K, So, Aro, and Aipo such that for any choice of r,r' 
there exists a projected (So, Aro, Aipo) -regular path from r to r' in T with < K 
collisions. 

Lemma 6. There exist M , So, Aro, and A(po such that given any (r, v±), (r' , v'j_) G 
C there exists a sample path from (r,v±) to (r',v'j_) with < M collisions and the 
corresponding projected path is (So, Aro, Aipo) -regular. 

The next two propositions deal with densities and are used in the proof of 
Prop. H 
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Proposition 4. Given Ar < Atq, Av± small enough, there exists t)n such that 
if \x is the uniform probability measure in the neighborhood (r — Ar, r + Ar) x 
(v± — Av±, v± + Av±) of (r, v±) £ C, then n > t)nv, where v is the uniform 
probability measure on C . 

Proposition 5. There exist r\, Ar, and Av± such that for any (r,v±_) £ C, 
T~'*o~(r.v ± ) > W{r,v±), where V( r ,v ± ) is the uniform probability measure onA( r v± ) = 
(r" - Ar, r" + Ar) x (v'[ - Av±,v'[ + Av±) for some (r" , v'[) £ C. 

Remark 3. In order to satisfy the conditions of Prop.[3]for both N th and (N+l) st 
push forward of <P, we first note that Prop.[3]holds similarly if (r, v±), (r 1 , v' ± ) £ C, 
where 6™ in = cos 3 (f - flo)^" and v™ ax = ^A^) since v ± in and V T X are 
arbitrary in the first place. Therefore to produce an N + 1 segment path from 
(r,v±) £ C to (r',v' ± ) £ C, one takes a one step regular sample path from 
(r,v±) £ C to some (f,v±) £ C and then follow an N step regular sample path 
to (r',v'j_). Prop. S] easily applies for N + 1 step regular paths. See proof in 
subsect. 15.41 

5.1. Proof of Lemma [31 Our periodic bounded horizon billiard configuration 
with the maximal flight distance denoted by T max has the following geometric 
property: 

Lemma 7. There exist Si, Ar\, and Atpi such that for each r £ UidDi there 
exists a 1-segment projected (Si, Ar\, Atpi)- regular path 7 that originates at r. 

Lemma [3 is a rather simple consequence of continuity and compactness. We 
include a proof for completeness purposes at the end of this subsection. 

Definition 4. Given 6%, Ar\, and Aipi from Lemma for any r £ dDi let 

U r be set of points in dDi that can be reached from r by a 1-segment projected 
(<5i, Ar\, Aipi) -regular path. 

Each U r is non-empty since r £ U r by Lemma [7] and open by definition of 
a projected (6%, An, Ziy?i)-regular path. Let J r be the connected component 
of U r that contains r. Then U r J r is a covering of a compact set dDi. Let J = 
{ Ji, • • • , J ni } be a finite sub-covering such that no element of J is fully contained 
in the union of some other elements of J . The later minimal property guarantees 
that each point in dDi is covered by at most 2 elements of J . Let r\, - ■ ■ ,r ni be 
the set of points that generated the original Ji :— J rj in the subcovering and 
rt j " " j be the set of midpoints of all the non-empty intersections of elements 
of J. 

Further re-order J±, ■ ■ ■ , J ni and f\, • ■ ■ , f ni clockwise such that fk £ Jk^Jk+i- 
Then there exists a projected (S\, Ar\, Z\</?i)-regular path 7$ : r\ — > si — > f\ — > 
Si f2 ■ ■ ■ — > r,, H — > s ni — > f ni — > s„ t — > r\ around the thermostat with 
Si, s\, ■ ■ ■ , s ni , s ni chosen to satisfy the (S±, Ar±, Z\(y9i)-regularity property. 

Assume now that r,r' are any two points on dDi. Then for some I and 
to, r,ri £ Ji and r',r m £ J m and there exist projected (S±, Ar\, Z\(/3i )-regular 
paths r — > s — > ri and r m — > s' — > r' . Therefore we can construct a projected 
(Si, Ar\, Zi^i)-regular path from r to r' that starts with r — > s — > ri, follows 7^ 
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from ri to r m , and ends with r m — > s' — > r' . The number of total collisions for 
such a path is bounded by Ki — 4n.j + 4. 

For any two thermostats there is a finite chain of length at most < p of 
thermostats joining these two such that any two neighbors in this chain can 
be connected by some 1-segment projected paths that do not meet thermostats 
tangentially. Let {7,} be a finite collection of such 1-segment projected paths 
such that the above connectivity property is satisfied using paths from 
only. Let 62 = ^ — max^H^I}, where the maximum is taken over all incoming 
and outgoing angles ip of all 7^. And let Ar2 and A(f2 be such that all 7, are 
(<5 2 , Ar 2 , Z\</?2)-regular. 

Then if we choose Sq = min{<5i, 82}, Aro = min{Z\ri, Ar 2 }, Aipo = min{Aipi, A1P2}, 
and K = J2i<i< P {Ki} where p is the number of thermostats in the system. 
Then there exists a projected (5q, Aro, Zl^?o)-regular path from any r € Ui-Dj to 
any r' Uj £>; with the number of collisions bounded by K . 



Proof (of Lemma^. Assume for S n = Ar^ — Atp^ = i there exists a 
sequence r n such that no 1-segment projected path originating at r n is (— , — , ^)- 
regular. Let r be an accumulation point of r n . Then by continuity there exists a 
1-segment projected path originating from r that is i ^)-regular for some 
m. Since (— , — , ^-regularity is an open condition, for each point f in a small 
neighborhood around r, there exists a 1-segment projected (-£-, -M-regular 
path originating from f , a contradiction. 



5.2. Proof of Lemma\d[ Suppose we start with r and move to r' along the pro- 
jected (So, Aro, Z\(/?o)-regular path from Lemma [5] Given v±, let oo be the cor- 
responding sample path that starts at (r, v±) and ends at (r',v±) for some vj_. 
We would like to extend <tq by o\ such that o~\ starts at (r',Uj_) and ends at 
(r',v' ± ). 

First of all, since the path in Lemma [5] is (So, Aro, Aipo)-regu\&i with < K 
segments v± has the following bounds: 

cos*(^ - 5 )if n < cos K £ - 6 )v ± <v ± < K ^ — < 



2 uy " " cos^(§ -<5 ) " cos*(f -S ) 

For simplicity, let us remove primes and tilde's for a moment and assume 
that 7 is a projected (So, Aro, ^<^o)-regular path from r back to r such that the 
corresponding sample path given v± maps (r, v±) to (r, vj_). The path 7 is easily 
produced by following several (So, Aro, Atpo)-regvlar segments forward and then 
the same segments backwards. It is expected that if we perturb 7 keeping the 
return to r, the final value of the normal velocity v± would deviate from the 
initial value v±. Repeating the perturbed path several times accumulates the 
deviation from v± with ultimate possibility to reach any values for the final 
normal velocity. 

Note that if 7 consists of only 2 segments the perturbation would not give v± 
different from vj_ since we have a restriction for the perturbed path to return to 
the initial point. Nevertheless 4-segment paths have enough flexibility to produce 
the desired effect. 
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Each 4-segment projected path 7 near 7 is uniquely determined by the initial 
position r and the four angles ipi, tf2, f3, and ^4 at which the particle reflects 
from the boundaries of the thermostats. For a path like that 

v± cos(y / 1 ) cos(y/ 2 ) cos(^) cos(y 4 ) 

Jrh) = JrVPl,<p2,<p3,<P4) = = 7 J / 1 7 x / s, 

V± COS((fl) C0S((y9 2 ) COS((/33j COS{<f4) 

where primes correspond to the angles of incidences for the corresponding 
segments. The constraint we have on / is that the path returns to the original 
point r. Note that / r (7) = 1. 

Lemma 8. For any r 6 dT there exists a r < I such that for any a £ (a r , 1], 
there exists a 4-segment projected (So, Aro, A(fo) -regular path from r to r with 
a property that for any corresponding sample path going from (r, v±) to (r, v±), 
we have ^ = a. Similar path exists with ^ = ~ property. 

Proof. First we are going to construct 7 such that f r (j) = 1 and there exists 7 
near 7 such that f r (j) 7^ 1 • 

Let 71 be a 1-segment projected (So, Aro, Z\(/?o)- r egular path originating at r 
and ending at some r' . Assume r € dDi and r' € dDj. There exists a segment 
72 that originates at r and hits Dj tangentially at some r|, possibly intersecting 
other thermostats. Since 71 is (So, Aro, Z\< / 9o)-regular, there exists r" £ dDi such 
that 7' := r — >• r' —> r" — > ri —> r is a 4-segment path having the first two 
segments (Jo, Aro, Ziy)o)-regular and not intersecting any other thermostats. If 
7 := r — > r' — > r" — » r' — » r, then / r (7) = 1 and f r (l') — 0, implying 
that / r restricted to paths r — > r' — > r" — > f — > r is an analytic function 
of one variable f £ Di and is not constant on the connected component of 7. 
Since (So, Aro, Z\</Jo)- r egularity is an open condition, there exists a projected 
(So, Aro, Z\</5o)- r egular path 7 := r — >• r' — >• f — >• r' — > r such that / r (7) 7^ L 

It may happen though that f r conditioned to return to r has a critical point 
at 7, in which case we are not guaranteed an open interval [a r ,l/a r ] of 
variation. However, since f r conditioned to return to r is not constant, there 
exists a projected (So, Aro, Aipo)-Tegulai path 7 near 7 such that f r conditioned 
to return to r does not have a critical point at 7. The conclusion of Lemma [5] 
follows. 

Given a < 1, let R a be the collection of initial positions r 6 Uj£)j such that 
for any a € (a, 1), there exist two 4-segment projected (So, Aro, Aipo)-~regu\a.r 
paths from r to r with -[JA. = a and = ^ respectively. Note that this time the 
4-segment paths are allowed to reflect off any thermostats. Each R a is open and, 
by Lemma[S]and continuity, {R a } ae ( ,i) i g a covering of the compact set UjDj. 
Extracting a finite sub-covering guarantees an upper bound a r < a max < 1 for 
all r £ Lii A- 

Therefore for every r £ UiDi and a £ (a ma x, ^ — ), there exist two 4-segment 
projected (So, Aro, Aipo)-reg\ileLT paths along which v± changes by a fraction a 
and A respectively. 

Let n be such that a max < [— — — ]" and be the if -segment 

v j_ 

projected (#0, ^?~0) ^</ J o)- r egular path from r to r' guaranteed by Lemma[5j The 
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corresponding sample path o~o given and v± maps (r, v±) to (r 1 , for some 
v±. Let 

a = mm] — , } " 

v± Vj_ 

Then a max < a < 1 since cos K (f - S Q )if n < v ± < co JJ^_ So) and v™ n < 
v '±_ < C™ ax . Let 75 be the 4-segment projected (So, Ar 0l Aip )-regul&r path from 
r' to r' such that — = 5 if v j_ < v± or — = i if v ± > w j_ . Then if we extend 
Ik by 75 repeated n times we obtain a projected (8$, Ar$, AipQ)-Yegv\ai path 
with the corresponding sample path a given v± going from (r, v±) to (r',v'j_) as 
desired. And the total number of collisions of such a path is < M = K + 4n. 



5.3. Proof of Proposition^ Lemmas O and [HI guarantee that there is a projected 
((5o, Aro, Z\</?o)- r egular path from (r, v±) to (r', with < M collisions provided 
that v*l in < v±,v' 1 _ < vf™. We would like to find an N > M such that for 
any choice of (r, v±) and (r',v' ± ) we can construct a path that has exactly N 
collisions. 

Lemma [7] provides us with a simple way to add two segments to a projected 
(<5o, Aro, ^v ? o) _r egular path with return to the same point r' with the same v'j_. 
However, there may be a parity issue if some of the paths from Lemma [S] have 
odd number of segments and others even. 

For each odd k, let Ck be the set of all r € UidDi such that there exists a 
projected (So, Aro, ^V 3 o)- r egular path from r to back to r that has k collisions. 
Each Ck is open. 

Claim. UfeCfc is an open covering of UidDi provided that So, Aro, and Aipo are 
small enough. 

Proof. First of all, we would like to choose So, Aro, and A(po, smaller than before 
if necessary, so that UfcCfc is non-empty. Let r € dDi and r' £ dDj be the closest 
two points connected by a projected path, i.e. \r' — r\ — r m ; n . Then there exists 
a rectangle with two of the sides parallel to the segment r — > r' that meets no 
other thermostats. Extend the sides perpendicular to r — > r' until they meet 
some other thermostats and Di at points ru and r; respectively. If r^ is at 
a corner of the rectangle, perturb r £ Di, r' £ Dj and/or ru £ -Dfc such that 
7~ — > f ' — > ffe — > r is a 3-segment projected path with all angles of incidence and 
refection not equal to ±-|; otherwise this property is guaranteed without the 
perturbation. Choosing (5o> Aro, and Aipo smaller than before if necessary, we 
ensure that this path is (So, Aro, Aipo)-regu\ai. 

Assume now that some r\ £ UjCj and let r^ £ Ck for some k. Then Lemma 
[5] provides a projected (So, Aro, Aipo)-Tegul&r path from n to r2 to which we 
append a path from r2 to itself with k thermostat collisions, and then the 
Lemma [5] path in reverse direction, i.e from r2 to r\. This generates a projected 
(So, Aro, Aipo)-Tegu\&r path from r% to itself with odd number of collisions with 
thermostats, a contradiction. 

Since U^-Cfc is an open covering of a compact set UidDi, there exists a hnitc 
sub-covering {Ck l: ••• , Ck m } with ki < fc 2 < ■ ■ ■ < k m . Therefore for each r £ 
UiDi there exits a fc m -segment projected (S , Aro, ^VcO-regular path 7™. 
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Let us choose N to be even. Then if the path from r to r' given by Lemma 
[5] contains odd number of segments, append it by 7™ to produce a projected 
(Sq, Aro, Zi</?o)-regular path with < K + k m segments. Then continue with boost- 
ing vj_ to the desired value v' ± as in Lemma |5] This produces a sample path a 
such that corresponding projected path has even number of segments and is 
(<5o, Aro, Aipo)-Tegu\&r with < N segments for some uniformly chosen even TV. If 
such a path happens to make < N thermostat collisions, we just add to it the 
required number of back and forth 2-segment projected (60, Aro, A(po)-regvlai 
paths guaranteed by Lemma [7] 



5.4. Pushing Density Forward: proof of Prop. [7} Given Ar < Aro and Av± < 
1 -mm cos N^ii_fi o ^ i e j. ^ k e un if orm probability measure in the neighborhood 

(7- — Ar, r + Ar) x (v± — Av±, v± + Av±) of (r,v±) £ C. We would like to show 
that 11 > f)NV, where v is the uniform probability measure on C. 

The projected (So, Ar , Z\<po)-regular path provided by Prop.[3]does not leave 
the set C max = {(r, v±) : cos N (% - S )v^ iD < v± < co JJT_ So) }- In the following 

Lemma [S] we will give the desired bound for any two points in (7 max that can 
be connected by a 1-segment (So, Aro, Aipo)-regu\a.r path. Prop. 0] follows by 
applying Lemma times. 

Lemma 9. Given Ar < Aro, Av± < ^v™ ln cos N (^ — Sq) and r\ > 0, there exist 

rf , Ar' , Av' ± such that if(r,v±), (r 1 ', Uj_) £ C max can be connected by a 1-segment 
(So, Aro, Aifo) -regular path and p is the uniform measure in the neighborhood 
(r — Ar, r + Ar) x (v± — Av±, v±_ + Av±) of (r, v±) £ C having density rj with 
respect to Lebesgue measure m on Q, then V*p > rf ft! , where p! is the restriction 
ofm to (r 1 - Ar',r' + Ar') x (v' ± - Av' ± ,v' ± + Av' ± ). 

Proof. When we push p forward by we first perturb it in the angle direction 
after which we obtain a product measure in 3-dimensional space with density 
T) x pv x (<p) at (f,<f,v±) £ (r - Ar,r + Ar) x (-|, |) x (v± - Av±,v± + Av±). 
If we only allow ip to vary by at most Aipo < then the lower bound on this 
density is r\ x p mm , where p m in is the minimum value of p$ ± ((p) when cos — 
S )vf n - Av±_ <v ± < co JJT_ 5o) + A Vl _} and \<p\ < § - S + A(p . Next we 
push this perturbed 3-dimensional measure forward by enhanced billiard map 
i.e. the billiard map with v± as an additional variable. The derivative matrix for 
such a map is 



/ TK+cos(y) r q 

cos(ip') cos(ip') 

TK,K fj rK cos(tp')-\-K f cos(t^) r At' +cos(tp' ) q 

cos (p' cos(^') 

, sin(y')^- ~ sm (y')^ cos(y)+cos(y') sin(y) cos(y') i 

\ COs(t^) COS 2 (ip) COS(^) / 



(8) 



and has determinant —1 (note the cancelation of the usual billiard map derivative 
CQh fr,\ with the perpendicular velocity change coe y) ). Therefore the lower bound 

COS(<£> ) L L J ° cos(^) ' 

on the density does not change when pushed forward. Now we need to project 
this pushed forward 3-dimensional measure back to (r, ujj-plane ensuring certain 
lower bound on the density of the projected measure. The only thing we need 
to guarantee is 
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Claim. There exist Ar' , Av' ± , Aip' that do not depend on (r,v±) and (r',v' ± ) 
such that the pre-image (under the enhanced billiard map) of the Ar' x Aip' x 
Z\w^_-neighborhood of (r' , <p>' , v' ± ) maps into Ar x Aip x Z\vj_ -neighborhood of 

(r,ip,v±). 

Then the lower bound on the density of V*\i in (r 1 — Ar, r' + Ar') x (v' ± — 
Av' ± ,v' ± + Av' ± ) is rf = r\ x p min x 2 Aip'. 

Proof (of Claim). Given Ar' < Ar^ and Aip' < Aipo, let A be the set of all 1- 
segment projected (backward) paths originating at f' at angle ip' with (r', <^') G 
(r' - Z\r', r' + Ar') x (<p' - A^', ip' + Aip'). We will denote the landing point of 
such a path by (r, ip). Also let 



G x {Ar' ,Ap>) = sup|f 



A 



G 2 {Ar',Aip') = sup\ip-ip\ 



A 



^ / a 1 a /\ 1 cos(<5) cos(os) 

G 3 (Ar',Aip') = sup 1 



,4 COs((/5') COs(y') 

G4( ^,^0=sup^- 

A COs(^') 

Note that for 1 < i < 3, G 4 (Z\r',ZV) uniformly as (Ar',Aip') ^ as 
long as all angles are bounded away from ^, which is ensured by the fact the 
r' — > r is a projected (#o, Z\r , Z\(/?o)-regular path. G4 is bounded by 1 < G4 < 
1/ cos(§ — (5 ). Therefore we can always choose Ar' and Aip' to ensure that 

Gi(Z\r',/V) < Z\r 



G 2 (Z\r', A/) < A/? 
G 3 (4r',A/K iax < ^ 

Now 

sup \v±- Vl _\ = sup I £2^^ £2^^ J < G 3 (Ar',Av')v'±+G4Ar',Av')Av' ± 

A A COS{ip') COS(ip') 



and we want it to be less than Av±. Choosing 

Av' ± < 



Av± 



2G 4 {Ar',Aip') 



does the job. 

Then for any f' E (r' - Ar' , r' + Ar'), p' E (tp' - Aip', ip' + Aip'), and v' ± E 
(v'j_ — Av'j_,v' ± + Av'j_), we guarantee that the image of (f',<p',v' ± ) under the 
enhanced billiard map belongs to the Ar x Av± x Aip- neighborhood of (r, <p>,v±). 
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5.5. Acquiring density: proof of Prop. [5j We would like to establish that there 
exist rj, Ar, and Av± such that for any (r, v±) G C, V^S^ r v± j > r\mA {r v±) , where 
m is Lebesgue measure on f2 and Au VJ _ \ = (r" — Ar, r" + Ar) x (v'\_ — Av± , v'[_ + 
Av±) for some (r" ,v'{_) G C. 

The 1-step push forward measure V*S/ r>vx \ is supported on a finite union of 
curves L in the phase space fl and is determined by the distribution p v± (tpi)dipi 
of the angle ip±. It may happen that (r, is a pre- image of (r', for several 
different values of <p\. Counting just one we conclude that the density at each 
(r',v'j_) is 

1 , \ 

> , === X PvAflh 

where J^j- and are entries (1,2) and (3,2) respectively of the enhanced 
billiard derivative matrix ©. 

The 2-step push forward measure 7*<J(r,uj_) is? m addition, determined by the 

distribution p v > ± (tf2) of </5 2 , where v' x = ^f$ v±. The density d(Vl8( r ^ ± )) of 
^*^(r,«±) at the endpoint (/',«") of the sample path starting from (r, «x) and 
determined by ipi and </?2 satisfies the estimate 

( dr" dr" 
dvj dv'l 
d<pi dip 2 

where the derivatives are computed using (O and a similar derivative matrix 
for the second flight. In particular, = ^p-§^ and = J^-- 

Lemma 10. There exists a uniform lower bound r\ on d(V^5t rtVx )) at any end- 
point (r",v'j_) of a 2-segment sample path for which the corresponding projected 
path is (So, Aro, A(po) -regular and initial v±_ is such that u™ ln < v± < u™ ax . 7n 
addition, the same bound holds at the corresponding endpoints if we allow the 
angles ipi and if 2 vary within Aip$. 

Note that at least one such path exists for any r by Lemma [7] 

Proof. By observing the entries of ((SJ we see that the determinant above is 
proportional to v±, which is bounded above by w™ ax . In addition, the terms in 
denominators are products of cos(</?i), cos(yi), cos(<^2), and/or cos(<p 2 ) of total 
power no greater than 3. Therefore by (Sq, Aro, Z\<po)- rc gularity of the path, the 
determinant is bounded above by some constant C — C(5q, w" lax ). 

In addition, if we let ipi and ip2 vary by at most Atpo, we get a similar 
bounding constant. The product p v± (ipi) x p v ' ± (if2) is clearly bounded for angle 
variations within A<po and i/™ m < v± < v™ ax . The conclusion follows. 

Unfortunately Aipo variation in <p\ and tf2 does not always produce nonzero 
variation in both Ar" and Avj_. For example, if we start from any location 
r G dr and choose r" — r, then we get no variation in v'[_ at r" . Note that 

det I e ^ avf ) =0 in this case. Indeed, such paths correspond to critical 
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points of v^_{ipi,tp2) constrained by r"(<^i,<^2) = const. In order to guaran- 
tee nonzero variation in r" and v'J_, for each r, we must find a 2-segment pro- 
jected (Sq, Aro, Z\(/?o)-regular path such that the above determinant is uniformly 
bounded away from zero. 

Lemma 11. For any r there exists a 2-segment projected (Sq, Aro, Aipo)-regular 
path 7 defined by angles of reflection ipi and if2 with | J(r, ip\, <f2)\ > 0, where 

J(r,tpi,ip 2 ) :=det U-H/^ 

Proof. Fix the first segment of 7 to be any 1-segment projected (So, Aro, Aipo)- 
regular path ending at some r' with angle of reflection ip\. Then J(r,tfi,ip2) is 
an analytic function of one variable (f2 ■ Similar to the proof of Lemma [5] there 
exist two 1-segment projected paths 71 and 72 originating at r' and ending at 
the same thermostat such that 71 is (Sq, Aro, ^A<A))-regular and 72 has the angle 
of incidence if = ±-|. Therefore, J(r,ipi,Lp 2 ) is not constant on the connected 
component and there exists ip 2 such that J(r, ipi, tp 2 ) > 0. 

For any r, let f(r) = sup{ J(r, tp 1 ,ip 2 )}, where sup is taken along all 2 segment 
projected (So, Aro, ^Vo)- regular paths. / is an everywhere positive continuous 
function on a compact set and therefore there exists the minimum value fo > 0. 
For each r, fix a 2-segment projected (So, Aro, Atpo)-Tegul&r path j(r) given by 

y>i(r) and^ 2 ( r ) such that | J(r,^ 1 (r),^ 2 ( r ))l > By the Lagrange multipliers 
method r"(ip\, tp 2 ) constrained by v'J_(ipi, (fi 2 )/v± = const does not have a critical 
point at (Tp 1 (r),Tp 2 (r)) and v'\_ (ipi, <p 2 )/v± constrained by r"(ipi,ip 2 ) — const 
does not have a critical point at (jp l (r), Tp 2 (r)). 

Let e(r) > be the maximal allowed variation such that any (f" ,v'[) £ 
(r" — e, r" + e) x (v'± — v±e, v'j_ + v±e) is an endpoint of a sample path given by 
(r, (p-L , if 2, v±) with ifx £ (Tp 1 (r) - Aip , ( r ) + Aipo (W 2 ( r ) - A Va , W 2 ( r ) + 

Aipo), and J(r, (fi(r), <fi2(r)) > e(r) is a continuous positive function on a 
compact domain dT and achieves its minimum somewhere. Let e = mine(r). 

r 

Then for each (r, v±) and r" and v'[_ given by j(r), at each (r",v'D £ A r v± = 
(r" - e, r" + e) x (v'[ - vf n e, v'{_ + vf a e), d(VJ [r , v±) ) > n. 



5.6. Proof of absolute continuity of the invariant measure with respect to Lebesgue 
measure m on fl. 

Lemma 12. // a is invariant under <P T and for all (r,v±) £ fl the absolutely 
continuous component ofV*St r ,vi_) * s nonzero, then u <C to. 

The assumption of Lemma [T^] follows from the proof of Prop. [5] for each 
(r,v±) £ Q, we constructed a projected (So, Aro, <4<£o)-regular path in some 
neighborhood of which the pushed forward measure is absolutely continuous 
with respect to m (we do not need a uniform lower bound on the density here) . 

Proof. Assume a is not absolutely continuous with respect to m. Then a can 
be decomposed as a sum of absolutely continuous and singular components u = 
/i<g + fi± with u^(f2) > 0. By the invariance of u 

M « + m± = » = V*hr,v±) = CP*"«)« + (7\V«)x + (?>*V)« + (P.VO-L (9) 
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Any push forward of absolutely continuous measure is absolutely continuous, 
which follows from the proof of Prop. |4]if we drop the regularity bounds and only 
keep track of absolute continuity. Therefore (P*H^.)±=0 and = CP*M<)«' I n 
addition, Prop.[S]guarantees that (V*fJ,±)± > 0, since we assumed that n±(f2) > 
0. Therefore 

^{Q) < + (7 > , a Mx)<(«) ) 

which is a contradiction since absolutely continuous part on the left hand side 
of ^ must be equal to the absolutely continuous part on the right hand side. 

5.7. Proof of Lemma QJ Let us first lift the density -^^v±e~ l3v2 - L to the 3- 

dimensional space by sampling ip from p v± (ip): the resulting density at (r, ip, v±) 
is 

2f3v ± e-^ Pv± (p). 

We push the measure with this density forward under the enhanced billiard 
map, the Jacobian of this map is equal to 1 (see Formula ©), so the value of 
the density of the pushed forward measure does not change at the corresponding 
new location (r' 3 tp',v' ± ). 

To find the value of the density of the pushed forward measure at (r', v' ± ) we 
just need to integrate over all possible values at (r, ip, v±) that map to (r' , ip',Vj_) 

as ip' ranges over (— -|, Note that v', — cos ( y ^ . Then we obtain 

J V 7T CO§ Z (<p) 

2" . OO 

2/3 J^J!il__ e ^K) 2 /co S 2 ( v ')^ = 2 ^ e ^W) 2 / J_ e -u* du 

V 7r cos z (p') J y'n 

-f ' -°° 

as desired. In the computation used a change of variables u — v^ w j_ tan(y>'). 
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